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By means of density functional theory (DFT) calculations (with and without inclusion of spin- 
orbit (SO) coupling) we present a detailed study of the electronic structure and corresponding 
microscopic Hamiltonian parameters of Na2lr03. In particular, we address the following aspects: 
(i) We investigate the role of the various structural distortions and show that the electronic structure 
of Na2lrC>3 is exceptionally sensitive to structural details, (ii) We discuss both limiting descriptions 
for Na2lrC>3; quasi- molecular orbitals (small SO limit, itinerant) versus relativistic orbitals (large 
SO limit, localized) and show that the description of Na2lr03 lies in an intermediate regime, (iii) 
We investigate whether the nearest neighbor Kitaev-Heisenberg model is sufficient to describe the 
electronic structure and magnetism in Na2lr03. In particular, we verify the recent suggestion of an 
antiferromagnetic Kitaev interaction and show that it is not consistent with actual or even plausible 
electronic parameters. Finally, (iv) we discuss correlation effects in Na2lr03. We conclude that 
while the Kitaev-Heisenberg Hamiltonian is the most general expression of the quadratic spin-spin 
interaction in the presence of spin-orbit coupling (neglecting single-site anisotropy), the itinerant 
character of the electrons in Na2lr03 makes other terms beyond this model (including, but not 
limited to 2nd and 3rd neighbor interactions) essential. 

PACS numbers: 75.10.-b,75.10.Jm,71.70.Ej,71.15.Mb 



I. INTRODUCTION 

The electronic and magnetic behavior of layered 5d 
transition metal oxides^ has been a subject of intensive 
discussion in the last years. Particularly exciting has 
been the suggestion by the authors of Ref. 2 that hexag- 
onal iridates such as Na2lr03 are a realization of the 
nearest neighbor Kitaev-Heisenberg (nnKH) model: 

H^=2KS7S] + JS i -S j (1) 

This proposal is based on the premise that spin-orbit 
(SO) coupling is the most important energy scale for the 
description of these systems so that Ir 5d t2 g orbitals are 
written in terms of j e s = 1/2 and j e ff = 3/2 relativistic 
orbitals, with the Kramers doublet j e g = 1/2 represented 
by the operator 5=1/2. The combination of Kitaev and 
Heisenberg terms leads to a complex phase diagram with 
various magnetic and spin-liquid phaseiP^. Obviously, 
some of these properties can only manifest themselves 
when the Kitaev term dominates or is at least compara- 
ble to the Heisenberg term. Also, other possible contri- 
butions, such as magnetic anisotropy, ring exchange or 
biquadratic exchange, to mention a few, may alter the 
phase diagram and the properties of the model consid- 
erably. Most importantly, while the Kitaev-Heisenberg 
expression is the most general fully-symmetric expres- 
sion for anisotropic pairwise magnetic interaction in the 
second order in spin in the presence of SO coupling (just 
as the Heisenberg exchange represents the same in the 
isotropic non-relativistic case), it is not necessarily short 
ranged in the presence of considerable itinerancy. 



So far, essentially all analyses of the nnKH model 
for Na2lr03 have been performed in the localized limit, 
where an assembly of weakly interacting relativistic 
atomic orbitals is assumed to be a good starting ap- 
proximation. On the other hand, first principles cal- 
culations suggest considerable derealization of electrons 
over individual Ir hexagons building quasi-molecular or- 
bitals (QMOsfH The associated "itinerant" energy scale 
(the band width) is fts 1.5 eV, to be compared to the 
single-site spin-orbit splitting scakP (3/2)A » 0.7 eV and 
the Hubbard and Hund's rule correlation energy scale of 
U — Jh ~ 0.5—1 eV. This makes the entire premise of the 
nnKH model questionable. At the same time, it has also 
been pointed outP^ that the nnKH model with the addi- 
tion of the 2nd and 3rd neighbors Heisenberg interaction 
is easier to reconcile with the experimental data. Such 
relatively long-range exchange interaction is another hall- 
mark of considerable itinerancy (here and below, when we 
speak of itinerancy, we imply mostly derealization over 
Ir 6 rings, but not necessarily over the entire crystal). 

In the present work we revisit and discuss the valid- 
ity of both limiting descriptions for Na2lr03; itinerant 
(QMO picture) versus localized (j e ff =1/2 Kramers dou- 
blet). To this end, we perform a thorough analysis of the 
electronic properties of Na2lr03 within non-relativistic 
and relativistic density functional theory (DFT) and de- 
rive, using projection on Wannier functions, the rele- 
vant hopping parameters and show that QMOs are nat- 
urally obtained as linear combinations of Ir ti g Wannier 
functions. We discuss the relation between the quasi- 
molecular orbital and the relativistic orbital, j e g, repre- 
sentations and show that the behavior of Na2lr03 lies in 



between a fully localized and fully itinerant description. 
Finally the parametrization of the electronic bandstruc- 
ture allows us to provide realistic estimates for the model 
parameters in the localized nnKH model. We thus inves- 
tigate whether we are close to a regime where the Kitaev 
interaction plays a decisive role or not. 

Quite unexpectedly, we find that Na2lrC>3 is an exam- 
ple of a material where minor details of the crystal struc- 
ture can dramatically affect the electronic structure, and 
simple guessing of the band structure parameters, or esti- 
mating them from simplified crystallographic models (so 
far all model calculations for this compound were uti- 
lizing one or the other approach) can be exceptionally 
misleading. In fact some of the models energetically dis- 
cussed in the community, while of undeniable theoretical 
appeal, are not even qualitatively close to the actual pa- 
rameter range in Na2lr03. 

While this particular compound is very intriguing and 
has been enjoying extraordinary popularity lately, we 
want to emphasize that this strong dependence of the 
electronic properties on details of the crystal structure is 
an important result, whose relevance goes beyond specif- 
ically Na2lr03 and is likely true for many other materials 
based on honeycomb transition-metal layers. 

The paper is organized as follows. In Section [n] we 
review the crystal structure and magnetic properties of 
Na2lrC>3. In Section |HT| we provide details of the DFT 
calculations and the projector method. In Section |IV| 
we present the results of the electronic structure analy- 
sis without inclusion of spin-orbit coupling and analyze 
the role of the structural distortions in Na2ir03. In Sec- 
tion [V] we investigate the role of spin-orbit coupling and 
discuss the relation between the QMOs and the relativis- 
tic orbitals (j e ff)- In this context, we discuss whether the 
existing experimental situation can distinguish between 
the DFT description (with the resulting itinerancy) and 
localized (j c g = 1/2) models. We proceed with an analy- 
sis of the single-site magnetic anisotropy in Na 2 Ir03 and 
find it to be relevant (pure j c ff = 1/2 states do not have 
any single-site anisotropy). In Section VI we provide ab 



imfio-derived estimates for the parameters appearing in 
the Kitaev and Heisenberg terms in Na 2 Ir0 3 and discuss 
the validity of the nnKH model by considering the exper- 
imentally observed magnetic order and attempts to ex- 
plain it from a local point of view. Finally in Section |VII| 
we present our conclusions. 



II. CRYSTAL STRUCTURE AND MAGNETIC 
PROPERTIES OF Na 2 Ir0 3 

Na2ir03 crystallizes in the monoclinic space group 
C2/m (No. 12^P (see Fig. [TJ and consists of Ir honey- 
comb layers (Fig. [I] (b)) stacked along the monoclinic c 
axis (Fig. [I] (a)) with an in-plane off-set along a. Na 
ions occupy both the interlayer positions and 1 /3 of the 
in-plane positions at the centers of Ir hexagons. This 
structure can be visualized as proceeding from Nalr02 
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FIG. 1: Crystal structure of Na2lr03. (a) Projection on the 
ac plane and (b) projection on the ab plane. 



with a Cdi2 structure with triangular Ir02 layers. In 
these layers 1/3 of the in-plane iridium atoms are sub- 
stituted by extra Na, i. e., its formula can be written 
as Na(Na!/3lr2/3)02, which, multiplied by 3/2, gives the 
usual formula of Na2lr03p'. 

An idealized crystal structure of this kind corresponds 
to having all nearest neighbor (NN) Ir-Ir and NN Ir-0 
distances equal and Ir-O-Ir angles of 90 degrees. The 
experimental structure of Na2lrC>3 departs from the ide- 
alized case and shows a few distortions: (i) orthorhombic 
distortion that introduces inequality among NN Ir-Ir dis- 
tances and among NN Ir-0 distances, (ii) IrOg octahedra 
rotations that place O atoms on the faces of a cube con- 
taining an Ir hexagon (see Fig. 2 of Ref. [5]) and (iii) 
trigonal distortion which is a compression of the IrOg oc- 
tahedra in the c-direction that induces a departure from 
90 degrees of the Ir-O-Ir angles. In Section |IV| we will 
discuss the effect of these distortions on the electronic 
structure of Na2lrC>3. 
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FIG. 2: Possible antiferromagnetic patterns in a honeycomb 
lattice 



As shown by transport, optical and high-energy spec- 
troscopy studiep™, Na 2 Ir0 3 is an insulator with an 
energy gap E g of 340 meV. Magnetic susceptibility 
measurements indicate a Curie- Weiss behavior at high 
temperatures with a Curie- Weiss temperature Qcw 

1 1 f j K and an effective Ir moment /i c ff = 1.82/xb- 
Na2lr03 orders antiferromagnetically below Tm — 15 K 
with an ordered magnetic moment fJ, or d ~ 0.2/ie- The 
fact that Tjv is much smaller than Qcw may be a sig- 
nature of frustration, but it may be also caused by the 
itinerancy of Ir 5d electrons^ as will be discussed in Sec- 
tion Ell 

The magnetic pattern observed experimentally corre- 
sponds to a zigzag ordering, in contrast to the prediction 
of a stripe order by the nnKH modeP (see Fig. [2| . Re- 
cently, Chaloupka et al^ argued that such a zigzag order- 
ing can be also obtained by the nnKH model, when one 
correctly includes all the terms contributing to NN Ir-Ir 
exchange. In Section |Vj we will discuss this proposition 
in more detail. 



III. METHOD 

In this work we perform DFT calculations using the 
linearized augmented plane wave (LAPW) method as 
implemented in the full-potential code WIEN2lf^. We 
employ the Perdcw-Burke-Ernzerhof generalized gradi- 
ent approximation^ to the DFT exchange-correlation 
functional and set the basis-size controlling parameter 
-R-KmaJ^ to 7. We consider a mesh of 500 k-points in 
the first Brillouin zone. Relativistic effects are treated 
within the second variational approach. Convergence 
with respect to relevant parameters (the k-point mesh, 
the RK max and the second variational energy cutoff, etc.) 
has been carefully checked. 



A. Calculation of hopping integrals 

In order to be able to discuss various Ir-Ir 5d processes, 
we parameterize our non-relativistic DFT results in terms 
of a tight-binding (TB) model where the TB Ir 5d hop- 



ping parameters are obtained through the Wannier func- 
tion projection formalism proposed in Ref. [15] and gen- 
eralized to molecular Wannier functions in Ref. [THl We 
first construct Wannier function projectors P^„(k) for 
the three i 2s Ir 5d orbitals and calculate the TB Hamil- 
tonian _ff TB (k) (in matrix form) via 

P TB (k) = P(k)#(k)pt(k), (2 ) 

where D(k) is a diagonal matrix of Ir 5d i 29 Bloch eigen- 
values and the matrix -P(k) is formed by the projectors 
P" lt/ (k). Here, indices a, m, and v run over equivalent Ir 
atoms in the unit cell, Ir i 2 g orbitals, and Bloch bands, 
respectively. Na 2 lr03 has two Ir per unit cell and only 
the six Ir i 2g bands near the Fermi level E-p are consid- 
ered in the construction of projectors. 

We calculate the hopping integral t™" l R Q /_ R / between 
orbital m on Ir atom a in the unit cell at a distance R 
from a reference unit cell and orbital m! on Ir atom a' 
in the unit cell at a distance R' from a reference unit 
cell by integrating H TB (k) over k- vectors in the first 
Brillouin zone: 

.mm' _ 1 rrTB /in -ik(R-R') (o^ 

"'a-R.a'-R.' — iy / , -"aro.a'ro' W K • \°) 

k k 

where ff™ Q , m ,(k) are the matrix elements of H TB (k). 
Correspondingly, the diagonal matrix elements t™™ give 
the on-site energies. 

B. Construction of quasi-molecular projectors 

As was argued in Ref. the most natural descrip- 
tion of the electronic structure of Na2lr03 is in terms of 
quasi-molecular (QMO) orbitals localized on a hexagon. 
The strongest Ir-Ir hopping is between 5d t2 g orbitals 
of neighboring iridium ions via common oxygens. In 
this case, an electron on a given Ir t2 g orbital propa- 
gates around an Ir 6 hexagon with the peculiarity than 
only a certain t2 g orbital at each Ir participates in the 
hopping^, e.g. Irl(xy)-Ir2(xz)-Ir3(j/z)-Ir4(a;j/)-Ir5(x2;)- 
lr6(yz) (see Fig. 2 of Ref.[5|. These QMOs are analogous 
to the molecular orbitals of the benzene molecule Ceile 
except for the fact that in benzene the same p-orbital 
on each carbon ion participates in the formation of the 
molecular orbital while in Na 2 Ir0 3 , as described above, 
different t 2 g orbitals are involved in one QMO and the 
three t 2 g orbitals on one Ir ion contribute to three differ- 
ent neighboring QMOs. We elaborate the details of the 
construction of the QMOs in what follows. 

QMO projectors P^ !/ (k) are obtained as linear com- 
binations of Ir i 2g projectors PMi/(k): 

P M Ak) = U MM T M (k)P Ml/ (k). (4) 

M 

where in the Ir £ 2g projectors Pjvf;,(k), the index M com- 
bines now the atomic index a and orbital index m, i.e. M 
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runs over all t 2g orbitals of all equivalent Ir atoms. With 
QMOs ordered as M — A lg ,E 2u , Elg, B lu , E lg , E 2u and 
Ir t 2g orbitals ordered as M = xy 1 , xz 1 , yz 1 , xy 2 , xz 2 , yz 2 
(the upper index labels Ir atoms), U is given by [u = 
exp(z7r/3)] 

/l 1 1 1 1 1 \ 

1 4 2 i 5 

1 LU LU —1 UJ LU 

1 LU 2 LU 4 1 U! 2 LU 4 , . 

^ = 11 1-1-1-1 ' (5 -' 

1 LU 4 LU 2 1 LU 4 LO 2 

\1 uj 2 lu 4 -1 lu 5 lu J 

and Tm(1c) are the Bloch factors, accounting for the fact 
that the 6 sites forming a QMO belong to several different 
unit cells. Actual values for these factors depend on the 
manner in which a particular band structure code selects 
the unit cell (see the Appendix for the WIEN2k settings) . 



Regarding hrst NN, we denote t\o and t\ q the hoppings 
between unlike t 2g orbitals via Op states (Fig. [3] (a)). t la 
and t\ a denote NN direct hoppings of cr-type. ijn and tjn 
denote NN hoppings between like orbitals lying in par- 
allel planes. In the ideal structure such hoppings consist 
of linear combinations with equal weight of ddir and ddS 
bonds. ti± and tjj_ denote NN hoppings between unlike 
orbitals lying in perpendicular planes (see Figs. [31(b) and 
(c))- 

Regarding the second NN hopping parameters, t 2 o 
and t 2 o denote hoppings between unlike orbitals via O p 
and Na s states (Fig. [3](e)). t 2a and t 2b (t 2a and t 2b ) de- 
note hoppings between like orbitals as shown in Fig. [3] (d) 
and t 2c , t 2 d and t 2e (t 2c , t 2( i and t 2e ) denote hoppings be- 
tween unlike orbitals (Fig. [3] (e)). 



A. Experimental crystal structure 



IV. NON-RELATIVISTIC ELECTRONIC 
STRUCTURE 

In this Section we analyze and discuss the Ir-Ir 5c? t 2g 
tight-binding parameters for Na2lrC>3 up to second near- 
est neighbors. As mentioned in Section [TTJ three struc- 
tural distortions are present in Na 2 Ir03: orthorhombic 
distortion, IrOg octahedra rotation and trigonal distor- 
tion. Besides, the stacking of the honeycomb planes 
inherently violates the rhombohedral symmetry even if 
each plane is ideal. The formation of QMOs relies on 
the dominance of mirahexagon hopping 5 and therefore 
is sensitive to structural details. Therefore it is impor- 
tant to understand the role of structural distortions in 
establishing electron hopping paths. This motivates us 
to study electronic properties of a number of artificially 
idealized Na2lr03 unit cells where structural distortions 
of different types are systematically eliminated^. Such a 
procedure has proven very usefuP^ in understanding the 
behavior of Sr2ix04. 

We consider four different crystal structures: (i) the ex- 
perimental crystal structure^, (ii) an artificially idealized 
Na2lr03 unit cell, Si, where the orthorhombic distortion 
has been removed from the experimental crystal struc- 
ture, (iii) an artificially idealized Na2ir03 unit cell, S 2 , 
where the IrOg octahedra rotations have been removed 
from Si, and (iv) an artificially idealized Na2lr03 unit 
cell, S3, where the trigonal distortion has been removed 
from 52. All TB hopping parameters calculated for the 
four structures are given in Table |T] and schematically 
represented in Fig. [3j We consider the following ratio- 
nale for labeling of the hopping parameters Eq. [3] In the 
experimental structure of Na2lr03 there are two first NN 
Ir-Ir distances and two second NN Ir-Ir distances due to 
the fact that the Ir 6 hexagons are not perfect. We denote 
the corresponding Ir t 2g - Ir t 2g hopping parameters as 
ti and t\ for the first NN and, respectively, t 2 and t 2 for 
the second NN hoppings. Further, we have various pos- 
sible hoppings between equal and different t 2g orbitals. 



Previous electronic structure calculations^ have iden- 
tified the dominant hopping integrals for Na 2 Ir03 to be 
tio and t 2 o [as well as ij and t 20 ; further on, if not 
explicitly stated otherwise, we refer to both equivalent ti 
(t 2 ) and i| (t 2 ) when writing ti (t 2 )}. In Table [I] column 
Sexp we present the complete list of hopping parameters 
up to the second nearest neighbors. A TB model based 
on these hopping integrals provides a good description of 
Na2lr03 Ir t 2g states near the Fermi level E-p [Fig. [4] (a)]. 

We first note a very good agreement between the ti o 
(~ 270 meV) and t 2 o (~ —75 meV) values obtained 
with our WIEN2k-based projection method and with the 
FPLO codeP as was used in Ref. [5] These leading Ir t 2g 
hoppings strongly tend to confine the electron's motion 
to a single Ir hexagon and, as a result, the electronic 
structure of Na2lr03 near the Fermi level is dominated 
by the formation of well separated and relatively weakly 
dispersive QMOs-. On an Ir hexagon, as shown above, 
each Ir atom participates with one of its t 2g orbitals (see 
Fig. 2 of Ref. [5j). These orbitals combine to form six 
QMOs according to the unitary transformation Eq. (J5J. 
In support of this picture, Fig. [5] (a) shows the density of 
states of Na2lr03 projected onto the six QMOs (singlets 
Ai g and Bi u and doublets E 2u and Ei g ), where states 
with certain predominant QMO character are clearly sep- 
arated in energy from one another. The near-degeneracy 
of Ai g and E 2u states around E-p is rather accidental re- 
sulting from the 0/^20 ~ —3.6 ratio (see Tablc|l]). The 
real-space representations of the QMO Wannier func- 
tions onto which the Na2lr03 DOS is being projected 
are shown in Fig. [6j The QMO Wannier functions were 
constructed as described in Section [TT] by explicitly ac- 
counting for the location of each Ir t 2g orbital in the 
crystal. 

Other NN and second NN hopping processes involving 
intraorbital and interorbital hoppings (see Table [I ) allow 
an electron to jump from one QMO to another and hence 
are responsible for the band dispersion. Many of those 
hoppings are of the same order of magnitude (although 
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FIG. 3: Schematic representation of Ir-Ir t2 g hopping paths up to second nearest neighbor in Na2lr03. 



mostly by at least an order of magnitude smaller) than 
<2 0, like, for example, t\\\ and tfu. For the "z" bond such 
hoppings will be between xz and xz or yz and yz orbitals 
(see Fig. [3](b)). These hoppings are equal to 47.7, 30.0, 
and 33.1 meV, depending on the NN bond (see Table]!]). 
In fact, such appreciable variations in magnitude, which 
violate the symmetry of an ideal Ir hexagon, are 
ubiquitous among the hoppings that connect neighboring 
QMOs. Some of them even change sign, as, for instance 
t\ a and t\ a . This feature results from the orthorhombic 
stacking, distortions within the b^Na planes, and rota- 
tions of IrOg octahedra. 



B. Structure Si obtained by removing the 
orthorhombic distortion 

We now consider an idealized Na2lrC>3 structure with- 
out the orthorhombic distortion of Ir hexagons; this 
structure, which we call S\ , as well as other structures in 
this Section, is tabulated in the Appendix. In the struc- 
ture Si: (i) all mtralayer Ir-Ir bonds are of the same 
length, i.e., the Dq^ symmetry of an Ir hexagon is re- 
stored, (ii) all NN Ir-0 bonds are of the same length, 
(iii) all Ir-O-Ir bond angles are equal to 98.7°, and (iv) 
the oxygens lie on the faces of a cube drawn around an Ir 



hexagon (see Fig. 2 of Ref. [5]). The 3D crystal structure, 
though, remains orthorhombic in this approximation, due 
to the presence of multiple Ir layers. This explains small 
residual variations among the nominally equivalent TB 
model parameters (Table[IJ column Si): E.g., comparing 
parameters labeled with and without overbar; also, onsite 
energies like the xy on-site energy is slightly lower than 
the xz/yz on-site energy. However, these variations of 
tig orbital on-site energies, as well as of equivalent hop- 
ping integrals, are now noticeably smaller than in the 
experimental Na2lr03 structure. 

We conclude that removal of the orthorhombic distor- 
tion restores (to a certain degree) the degeneracy of the 
Ir tig orbitals, but does not change the hierarchy of hop- 
ping integrals. In the structure Si, the t\o and tio 
values are close to the respective values in the exper- 
imental Na2lrC>3 structure and, as a consequence, the 
overall structure of the ti g bands is only slightly changed 
[Figs. |4| (b) and[5](b)]. 

C. Structure S2 obtained by removing the IrOs 
octahedra rotations 

In the structure Si that we designed in the previous 
Section, two types of distortions are still present: (i) trig- 
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a) Experimental structure 
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(b) Idealized structure S. 
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FIG. 4: Na2lr03 bandstructures near the Fermi level Ef = 
calculated using DFT (black solid lines) and a TB model that 
considers up to NNN hopping processes between Ir t2 g orbitals 
(dashed lines). The data are obtained with (a) experimental 
crystal structure and idealized structures (b) Si, (c) 52, and 
(d) S 3 . 



TABLE I: Nearest neighbor (NN) and second NN hopping in- 
tegrals between Ir tig orbitals for the experimental structure 
and three idealized structures Si, S2, S3 of Na2lr03 (see text 
and Appendix for a description of the structures and param- 
eter labeling). The NN = data are Ir ti g on-site energies 
and mierorbital hoppings; the NN = 1 and NN = I (NN = 2 
and NN = 2) data are hoppings over nonequivalent (due to 
orthorhombic distortion) NN (second NN) Ir bonds. 



NN 

~0~ 



xy — > xy 
xz — > xz 
yz — > yz 
xy — > xz, xy — > yz 
xz — > yz 

xy -> xy (ti||) 
xy -> xz, xy -¥ yz (t 10 ) 
xy -¥ xz, xy -> yz (tix ) 
xz — > xz, yz — > yz (ti\\) 
xz — > xz, yz — > yz {tier) 
xz — > yz (til.) 

xy -> xy (ti a ) 
xy -> xz, xy -> yz (t I± ) 
xz — > xz, yz — > yz (ti\\) 

xz -> yz (tio) 



xy - 
xy 
xy - 
xy 

xz - 
xz 



xy — > xy (t 2o ) 
» xz, xy — ¥ yz 
-¥ xz, xy 
» xz, xy 
-> xz, xy 
> xz, yz 
-» xz, yz 

xz 

xz 



> yz 
yz 

> yz 

yz 

> yz 



(ho) 

(fee) 

(M 

(he) 

(*2a) 



yz (t2e) 

yz (fed) 



xy -> xy (t,2 b ) 
xy xz, xy -»■ yz (t 3e ) 
xy xz, xy -»■ (t 2(i ) 
-> xz, yz -> 2/2 (t So ) 
-> yz (t 5o ) 
xz -¥ yz (t 5c ) 



-448.8 
-421.5 
-421.5 
-27.8 
-23.1 

47.7 
269.6 
-25.6 

30.0 
-20.7 
-21.4 

25.4 
-11.9 

33.1 
264.4 

-3.5 
-75.8 
-36.5 

12.5 
-21.4 

-0.6 

-1.5 
-18.6 

10.2 

-1.4 
-19.0 
9.3 

-1.4 
-77.0 
-30.4 



Si 
-422.9 
-421.8 
-421.8 
-26.4 
-25.2 

34.1 
268.5 
-16.6 

33.2 
3.5 
-16.4 

0.2 
-17.6 
33.9 
264.8 

-2.6 
-77.4 
-35.3 

10.1 
-19.2 

-3.1 

-1.6 
-19.0 

10.2 

-1.4 
-19.2 
10.2 
-3.0 
-78.0 
-35.1 



S 2 
-422.8 
-421.2 
-421.2 
-21.2 
-18.8 

27.8 
231.7 

43.7 

17.2 
-66.5 

41.7 

-65.5 
46.9 
21.2 

228.7 

-18.9 
-94.7 
-52.1 
1.7 
-7.3 
-16.6 
-1.0 
-7.1 
2.4 

-1.2 
-8.4 
0.7 
-17.7 
-95.2 
-51.6 



S3 
-601.1 
-601.1 
-601.1 
-13.5 
-14.7 

120.8 
209.7 

-5.3 
118.9 
-381.6 

-4.9 

-382.8 
-5.3 
120.5 
211.7 

2.0 
-82.1 
-38.5 
6.9 
1.9 
1.4 
5.7 
2.4 



5.7 
2.1 
7.5 
1.5 
-81.9 
-38.9 



onal squeezing of Ir06 octahedra along the (111) direc- 
tion perpendicular to Ir hexagon planes and (ii) IrOg 
octahedra rotations that place O atoms on the cube's 
faces. We now consider structure S2, where the IrOg oc- 
tahedra rotations are removed from S±. In this structure, 
the Na-0 and Ir-0 bond lengths are the same (in the ex- 
perimental structure, the former is considerably longer). 
This feature enhances the second NN hopping processes 
through Na s states, such as £20, t%ai tic (and the equiv- 
alent overbar hoppings) as shown in Table |lj column S2. 
At the same time, the NN O-assisted hopping txQ gets 
reduced and the t\ o/foo ratio decreases to ~ —2.4, re- 
sulting in a larger separation of the lowest (B\ u ) band 
from the rest of ti g bands [Fig. E] (c)]. Formation of 
QMOs still takes place in structure S2 [Fig-[5]( C )L but the 
QMO bands are more dispersive compared to the exper- 
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FIG. 5: Na 2 Ir0 3 DOS projected onto QMOs for (a) experi- 
mental crystal structure and idealized structures (b) Si, (c) 
52, and (d) 53. The Fermi level is set to zero. 



imental or S± structures, due to increased inierhexagon 
NN hopping integrals t\ a , t\\_ (and equivalent t\ a , ti_i_): 
thus, one observes broadening of the A\ g band and re- 
distribution of weight away from the Ei u doublet. 



D. Structure S3 obtained by removing the trigonal 
distortion 

We finally consider a most idealized Na2lrC>3 structure 
S3 without the trigonal distortion, i.e., with 90° Ir-O-Ir 
bond angles. Importantly, one can only remove this dis- 
tortion, while keeping the Ir-0 bond length the same, if 
the Ir-Ir bonds are shortened. Because of that, the hier- 




FIG. 6: Real-space representation of the QMOs in Na2ix03 
obtained by the Wannier projector method. 



archy of hopping integrals changes drastically (Table |TJ 
column 5*3). The dominant hopping is now the direct 
NN hopping between like orbitals t\ a (and the equiv- 
alent ti a ) reaching ~ —380 meV, while the O-assisted 
hopping tiQ (<io) nas been reduced to ~ 210 meV. Ac- 
cordingly, the large mterhexagon interaction destroys the 
QMO picture, as illustrated by the strongly dispersive ti g 
manifold in Fig. [4] (d) and the derealization of individual 
QMO characters over the whole DOS range in Fig.[5](d). 
We also observe that the main reason for the trigonal 
squeeze is the geometrical effect of optimizing simultane- 
ously the Ir-Ir and Ir-0 bonds. As a result, even though 
the on-site t2 g orbitals split into an a\ g singlet and an e g 
doublet, this is not a strong effect and not the driving 



force for the squeeze, as it is often assumed in the spirit 
of localized limit and the Jahn- Teller effect. 

Summarizing these results, in the S3 structure, the NN 
direct hopping increases by an order of magnitude com- 
pared to the experimental Na2lr03 structure and the NN 
O-assisted hoppings get suppressed. Therefore we con- 
clude that structural distortions of all types in Na2ir03 
act constructively to enhance the mirahexagon effective 
hopping parameters (such as t\ o and £20) an( i suppress 
the mierhexagon ones (such as NN direct hopping) fa- 
voring the formation of QMOs. 

V. SPIN-ORBIT COUPLING 



0.5 




Tt r x m r 



We proceed now with the analysis of the electronic 
structure of Na2lr03 in the presence of spin-orbit 
(SO) coupling. Previous relativistic DFT calculations^ 
showed that Na2ir03 states near the Fermi level expe- 
rience strong relativistic splitting with pronounced con- 
centration of jeff = I character in the upper two bands. 
However, the Na2ir03 relativistic states seem to preserve 
their QMO identity as well [see Fig. S6 (b) of Ref. [5]. In 
order to understand such duality, we set up a TB model 
for the Ir t 2g orbitals that includes also local SO inter- 
action terms. With this TB+SO model, we are able not 
only to confirm the relativistic DFT results by calculat- 
ing DOS but also to access the composition of individual 
states and trace their evolution as a function of the spin- 
orbit coupling A. 



FIG. 7: WIEN2k relativistic bandstructure (black solid lines) 
versus TB+SO model relativistic bandstructure (red dashed 
lines) of Na2lr03 as described in the text. In model calcula- 
tions, A = 0.44 eV was used. 



Since our purpose is to reconcile the QMO and rela- 
tivistic orbital (RO) pictures, we analyze the AL • S ma- 
trix elements between spin-f and spin- .J, QMOs to see 
how SO coupling mixes QMO characters. They can be 
easily obtained by applying the unitary transformation 
UT(k) [Eq. pH) ] to the AL • S matrix in the t 2g basis: 



-^qmo s ) 



UT(k)H? 2 °T H (k)U H . 



(7) 



This equation explicitly illustrates how k-vector depen- 
dence enters the SO matrix elements in the QMO ba- 
sis. Concise expressions can be derived if one notes 
that QMOs can be represented by their "winding num- 
ber" n which defines a phase change = ^ of t 2g 
orbitals around a hexagon. In this notation, QMOs 
non-relativistic DFT Ir t 2q bands of Na 2 Ir0 3 . It includes A ig > E 2u, E lgi B lu , E lg ,E 2u correspond to, respectively, 



TB+SO model 



We start with a TB model that perfectly describes the 



three hundred and twenty one hopping integrals between 
up to 50 nearest neighbors. We then double the dimen- 
sion of the TB Hamiltonian matrix to introduce spin de- 
pendence and add local SO coupling terms (AL • S) that 
mix spin-f and spin- 1 subspaces: 



n = 0,1,2,3,4,5 winding numbers. The AL • S matrix 
elements in the QMO basis are given then by 

o n A. in 1 -.).. M — n)ir , 
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Importantly, even though SO coupling is a local on-site 
interaction, it couples neighboring quasi-molecular or- 
bitals and therefore is k-vector dependent in the QMO 
basis. 

Having thus set up the TB model, we vary the SO cou- 
pling strength A until the best matching with the DFT 
relativistic bands is achieved, which is found to corre- 
spond to A = 0.44 eV (Fig. [7). 
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and 
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We list numerical values of the matrix elements for two 
representative k-vectors: k = (0, 0, 0) (point T) (first two 
tables) and k = (f-,0, 0) (last two tables). 
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Several comments are in place here, 
coupling mixes QMOs at all k-vectors 
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First, spin-orbit 
Even at the T 
point, i. e., on the same hexagon, the three upper QMOs 
{A\g and two i? 2 u) are SO coupled to the three lower 
QMOs (Bxu and two Ex g ), which explains sizable shifts 
of the relativistic bands compared to the non-relativistic 
ones at this k- vector. Additionally, SO coupling induces 
splitting of the degenerate Ei u and Ex g states at all k- 
vectors. Another striking feature of the calculated AL • S 
matrix is that its A lg , E 2u (upper triplet) and B lu , Ei g 
(lower triplet) blocks are identical. This means that if not 



for the accidental near-degeneracy of the A\ g and Em 
states (which magnifies the SO induced energy shifts) 
the upper and the lower triplets would have been equally 
affected by the SO coupling. 



B. Quasimolecular orbital basis versus relativistic 
basis 



The main difficulty in describing the Na 2 Ir03 band- 
structure is that it interpolates between eigenstates of 
two Hamiltonians: the itinerant TB Hamiltonian of (pri- 
marily) mtrahexagon electron hopping that preserves the 
s z spin subspace and the local spin-orbit (SO) interac- 
tion AL ■ S Hamiltonian that couples different spin sub- 
spaces. The eigenstates of the TB Hamiltonian are quasi- 
molecular orbitals (QMOs), while the eigenstates of the 
SO interaction (in the i 2g subspace) are relativistic or- 
bitals (ROs) |jeff,j|fj) characterized by an effective total 
angular momentum j' e g and its ^-projection j^: 
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This basis^ can be explained as follows; three t 2ff or- 
bitals (total degeneracy, including spins, is 6) are split 
into a lower-lying quartet j c g = 3/2 and an upper ly- 
ing jeff =1/2 doublet, and the 5d-electrons of Ir 4+ fully 
occupy the lower quartet leaving the upper j cf f =1/2 
doublet half-filled. This makes this situation similar to 
a non-degenerate Hubbard model (S=l/2 doublet on a 
site), with the important difference that in the Hubbard 
model the hopping matrix elements preserve the s z spin 
subspace, while here the states of the j e ff =1/2 doublet 
are spin-orbit mixed states, leading to a strong anisotropy 
of hoppings and their dependence on spin (or rather to- 
tal moment) direction. This may bring about anisotropic 
exchange, e. g., the Kitaev exchange on a honeycomb lat- 
ticP. 

By gradually increasing an effective spin-orbit coupling 
strength A c ff, 

A 2 

Aoff = r. „ , ,0 , h o = 0.270 eV, (15) 



(t 



1 o 



A 2 ' 



from to 1, one can trace a smooth evolution of the 
TB+SO model eigenvalues from, respectively, the non- 
relativistic (QMO) limit to the fully relativistic (RO) 
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FIG. 8: Properties of the TB+SO model of Na 2 Ir0 3 at the 
r point as a function of effective SO coupling A e ff defined in 
Eq. | |15| . The vertical dotted line marks the realistic A c fr = 
0.73 value for Na 2 Ir0 3 . (a) The eigenvalues of the TB+SO 
model. At all A e ff points except A e fi = 0, a A c fi dependent 
scaling has been applied to the eigenenergies to keep them 
within the [—1.5,0.2] eV range, (b) The j off = § (solid line) 
and Biu (dashed line) weights on the lowest state, (c) The 
jeff = 1 (solid line) and total Ei u (dashed line) weights on the 
uppermost state. Inset shows individual contributions from 
the two E 2u QMOs. 



limit (see Fig. [8] (a) for the data at the T point). An 
SO coupling parameter of A = 0.44 eV for Na2lr03 cor- 
responds to A c ff = 0.73, which is marked by a vertical 
dotted line in Fig. [8] 

The RO basis is an attractive starting point to describe 
the low-energy physics of Na2lrC>3 as it allows to trun- 
cate the Hamiltonian to only j e fi = | states that domi- 
nate near the Fermi energy and map Na2lr03 onto the 
Kitaev-Heisenberg model. Although this approach might 
seem reasonable given the noticeable separation of the 
jeff = \ and jeff = | characters in the DOS of Na2lr03 
[c/. Fig. 2 (b) of Ref. 122] . we argue that the itinerant 
terms are too strong to be neglected (which should not 



be surprising since A = 0.44 eV < W ~ 4iio = 1 eV) 
and that, consequently, the QMO basis is as well (or as 
poorly) justified to work with as it is the RO basis. 

To support this statement, let us concentrate on the 
TB+SO model states at the T point. Fig. [8] (a) shows 
the evolution of the model eigenvalues as a function 
of A G ff (Eq. 15 1. In the non-relativistic limit (A c ff = 
0), the states are almost purely (with slight devia- 
tion due to orthorhombic distortion) QMOs, ordered as 



Biu, Ei g , Ai g ,E2u with increasing energy^. At the same 
time, at each state the j c s — \ contribution is 1/3 and 
the joff = | contribution is, correspondingly, 2/3 (for 
one of the two Ir atoms). Note that, since the model dis- 
tinguishes spin-f and spin- J. states each level is doubly 
degenerate. 

With the QMO splitting obviously prevailing for zero 
SO coupling, we now want to quantify the QMO char- 
acter rectification upon increasing A c ft by calculating the 
QMO and RO weights on two selected states: the low- 
est (B lu ) and the uppermost (E 2u ). The B lu state 
[Fig.[8](b)] is a simpler case as it is non-degenerate (apart 
from spin) and quite well separated from the rest of the 
QMOs so that the SO effects here should be less impor- 
tant. Changing X e g from to 0.73 (Na2lr03 value), the 
Jeff = | weight on this state increases from 0.6667 to 
0.8320, whereas the B\ u weight is only slightly reduced 
from 0.9932 to 0.9567. This indicates that the lowest 
relativistic state at the T point in Na2lr03 is better de- 
scribed by a QMO Bi u than by one of the j e ff = § ROs. 
In fact, this turns out to hold for the whole lowest rela- 
tivistic band [cf. the j c ff- and QMO-projected Na2lr03 
DOS in, respectively, Fig. 2 (b) of Ref.|2Jand Fig. S6 (b) 
of Ref. [5]. 

The uppermost state is one of the E 2u doublet states. 
It is near-degenerate with A\ g and the other E 2u and, 
therefore, the SO effects are here particularly strong. At 
the r point, though, it can only couple to itself or to the 



linear combination of t 



other E 2u [see Eqs. (10) and (111], depending on which 



iese degenerate states is consid- 



ered. Upon switching A c ft on, the jeS, — § weight on 
this upper states rapidly grows from 0.3333 to ~ 0.6 in 
the range < A e ff < 0.05, and then gradually increases 
to 0.8295 at A cff = 0.73 [Fig. [8] (c)]. At the same time, 
the weight of one of the E 2u states (we may call it E 2u ) 
is reduced from 1.0 to 0.53730 [see inset of Fig. [8] (c)]. 
However, the total weight of two E 2u states is barely 
changed: at A ff = 0.73 it equals 0.9617. This means 
that the uppermost relativistic state at the T point in 
Na2lr03 is very well described by a linear combination 
of two E 2u states (which is also a QMO) with, in general, 
Aeff-dependent individual contributions. 

The Bi u and E 2u states (at A = 0) seem to simultane- 
ously bear both RO and QMO features up to very strong 
SO coupling, with the QMO character even dominating 
at Acff < 0.9. This can also be illustrated by inspecting 
the composition of, e. g., the lowest energy band state as 
shown in Table [TT] At zero SO coupling, the doubly de- 
generate lowest state corresponds to (almost) pure B\ u "f 
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and B lu | QMOs 25 . At A cff = 0.73, the structure of this 
state is strikingly similar to the B\ u states, with only 
slight admixtures of the xz and yz orbitals of opposite 
spin. Even as some very high A e ff, when the RO j e g — | 
weight is close to 1, the states retain the B lu f and B\ u \ 
QMO features. 

The features shown in this Section, not unexpectedly, 
characterize Na2lrC>3 as intermediate between the non- 
relativistic (pure quasi-molecular orbital) and fully rela- 
tivistic (pure RO) cases. 

Moreover, these results show that, in the RO repre- 
sentation, the upper band states are not pure j e g =1/2 
states but there is some significant mixing of j Q ff = 3/2 
states. In fact, for the upper band states, the projec- 
tions onto jeg = 1/2 and j off = 3/2 are, respectively, 
0.64 and 0.21 with 2(0.64 2 + 2 x 0.21 2 ) = 1, while in 
the non-relativistic case these projections are both equal 
to \ZT/6 = 0.41. Note that looking at the weights 
may be misleading. Indeed this state appears to be 
2x0.64 2 = 82% pure j eS = 1/2 state [Fig. |] (bottom)], 
but its projection on the j e ff = 3/2 state is only twice 
smaller than in the non-relativistic case. In other words, 
the hopping between the upper Kramers doublets, ini- 
tially not considered in Ref. is only reduced by about 
a factor of two compared to the non-relativistic case. One 
but possibly not the only consequence of this fact is that 
the contribution of the Kitaev term in the analysis below 
may be overestimated, probably by as much as a factor 
of two. 



C. Comparison with experiment: branching ratio 

An argument frequently used to justify the assump- 
tion of pure ROs in Na2lr03 is that it is experimentally 
supported. However, the experimental evidence is incon- 
clusive. It is first assumed that the electronic states are 
pure ROs and then it is shown that this assumption does 
not contradict the experiment, yet the experiments, upon 
a closer look, do not falsify the DFT picture, either. A 
typical and, by far, the most often used quantity to dis- 
cuss the nature of the states in iridates is the branching 
ratio (BR) extracted from X-ray absorption spectroscopy 
(XAS) experiments. In XAS, essentially, (L • S) is mea- 
sured. This expectation value is of course zero without 
spin-orbit coupling. A detailed and very insightful anal- 
ysis can be found, for instance, in Refs. in par- 
ticular, it is shown that, for a related iridate, the main 
contribution to (L • S) (1.4 out of 2.1) doesn't come from 
the t2 9 orbitals, which define the j c a = 1/2 states, but 
from the admixture of the e g orbitals. In our calculations 
-shown below- we observe the same behavior. 

We apply our TB+SO model to calculate (L • S) for 
Na2lr03 where L and S are, respectively, the total orbital 
and spin angular momenta of Ir 5d electrons. (L • S) is 
related to the experimentally accessible branching ratio 



BR = 



(2-Q 
(1 + r)' 



(L-S) 

n h 



(16) 



with nh = 5 being the average number of 5d Ir hole d 30 * 3 U 
In recent XAS measurements^, BR = 5.5 — 5.7, trans- 
lating to (L-S) = —2.7ft 2 , was obtained for Na 2 Ir03 and 
interpreted as a sign of strong spin-orbit coupling. 

When applying the TB+SO model that we constructed 
for Na2lr03 in Section VA the calculated (L-S) = 
— 0.73ft 2 (as compared to — lft 2 in the limit A G ff = 1). 
This value is several times smaller than the experimental 
value. This is, however, not unexpected given the signif- 
icant contribution of the Ir e g empty states to (L • S) (c/. 



Ref. [26]) . which are not considered in the TB+SO model 
discussed in the previous Section. In order to make a 
meaningful comparison with experiment, we extend our 
TB+SO model to include (in the same spirit) also the 
Ir e g states. (L • S) within such a model is — 1.91ft 2 . This 
is about 30% less than the experimental value reported by 
Clancy et alP 2 . This result is indeed in good agreement 
with experiment, given the large fluctuations in experi- 
mental values. For instance, Ref. [3H reported (L-S) = 
-3.1ft 2 for Sr 2 Ir0 4 while Ref. [Ml reported -2.1ft 2 (about 
30% difference) for the same compound. This example 
gives a sense of possible fluctuations between results of 
different experimental groups, and therefore our theoret- 
ical (L • S) value for Na2lr03 might be even closer to the 
true result. 

The main conclusion from these calculations is that 
with the TB+SO model based on all five Ir 5d orbitals we 
are able to reasonably reproduce the large experimentally 
measured (L • S) value in Na 2 Ir03, which validates our 
approach. As our analysis shows, the large (L • S) does 
not necessarily mean an ideal separation of j e ff = § and 
jeS — g RO states, but rather the effect of e g states also 
contributing in the process. Due to the peculiar electron 
hopping hierarchy in Na 2 Ir03, QMOs might be a better 
basis. 

In conclusion, the XAS experiments only tell us that 
the upper Kramers doublet has a considerable contribu- 
tion coming from j c g = 1/2, but not that it is a pure RO 
state. 



D. Comparison with experiment: RIXS 

Another experiment sometimes quoted as supporting 
the fully relativistic j D s = \ picture is resonant inelastic 
x-ray scattering (RIXSJP^ In this experiment a joint den- 
sity of electronic states (JDOS) is probed, somewhat sim- 
ilar to that in the infrared absorption but with different 
matrix elements. The authors of Ref. [33j observed several 
peaks in JDOS, of which the lowest peak at ~ 0.42 eV 
was interpreted as transitions across the Mott-Hubbard 
gap, consistent with a 30% smaller optical absorption 
threshold. The next two peaks are close to each other at 
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TABLE II: Expansion coefficients of the lowest doubly degenerate energy states of the TB+SO model in the ti g basis (The 
upper index of the t2 g orbitals labels Ir atoms in the unit cell). The coefficients are given for three A (A e fi) values. The Bi u 
and Jeff = | weights of the various states are given at the bottom of the table. 
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0.72 and 0.83 eV and were ascribed to transitions from 
the j e f[ = 3/2 quartet into the upper j c g = 1/2 doublet. 
The splitting of 110 meV was ascribed to the trigonal 
splitting. Altogether, this interpretation suggests an SO 
coupling A - §( 0-72+0.83 _M2 } _ 0.39, a very reasonable 
number, if slightly too small. 

This analysis, even though it looks reasonable on the 
first glance, has serious shortcomings. First, the deduced 
trigonal splitting is nearly twice as large as the actual 
trigonal splitting. In fact, the trigonal splitting is de- 
cided by the electrostatic field of the ligands, and by one- 
electron hoppings thereto; both are very well accounted 
for by the DFT calculations, which give At = 75 meV. 
Second, even a At = 110 meV cannot produce well sep- 
arated peaks in JDOS, given that the Ir-Ir hopping is 
tio = 270 meV. Third, even if one completely neglects 
the Ir-Ir hoppings^l, in order to extract A and At one 
has to diagonalize the full Hamiltonian including both 
factors and then fit the resulting eigenvalues to the ob- 
served peaks. After doing that, one gets A = 0.5 eV and 
At = 180 meV. The latter number is more than twice 
the actual trigonal splitting. 

On the other hand, this experiment is consistent with 
DFT band structure. To demonstrate that, we have per- 
formed DFT calculations for the magnetic zigzag phase. 
We note that the results do not depend qualitatively on 
the choice of the pattern and the magnetization direction. 
In order to account for the missing correlation effects 
and adjust the direct gap to be consistent with infrared 
measurements^, we applied a rigid shift of 200 meV 
between the occupied and empty bands ("scissor oper- 
ator"). This exercise gives a JDOS which has a broad 
feature, consisting of (i) a peak at 0.42 and a shoulder 
0.48 eV (compared to 0.42 eV in the experiment) corre- 
sponding to the transition between the top QMOs and 
(ii) a peak at 0.77 eV and a shoulder at 0.81 eV cor- 



responding to transitions from the lower QMOs. While 
the experiment finds two peaks at 0.72 and 0.83 eV, one 
should keep in mind that the matrix elements, omitted in 
our calculation, can easily suppress or enhance a shoul- 
der, making it disappear (at 0.48 eV) or become a sep- 
arate peak (at 0.81 eV). Therefore we conclude that the 
agreement between experiment and our calculations, sim- 
plified as they are, is reasonably good. 



VI. MAGNETISM 

A. Nearest neighbor Kitaev-Heisenberg model 

One term neglected in the conventional Kitaev- 
Heisenberg model treatment is the single-site magne- 
tocrystalline anisotropy. Localized electrons with the 
spin 1/2 do not have any anisotropy, no matter how 
strong the spin-orbit coupling is. However, if hopping 
is considered, electrons can have a preferred spin direc- 
tion, which in the language of the nnKH Hamiltonian 
would be reflected in a single-site term proportional, in 
the lowest order, to A • S where A is a vector. Such 
terms are usually neglected when dealing with the nnKH 
model. Our calculations^ without including U show a 
magnetic anisotropy as large as 3 meV per Ir (in order 
to address the single-site anisotropy, we compared ferro- 
magnetic calculations) . This energy should be compared 
to the total magnetic stabilization energy (i.e. the energy 
difference between magnetic and non-magnetic solutions) 
of maximally 5 meV. When the DFT calculations are per- 
formed including a U = 2 eV, the magnetic anisotropy is 
as large as 8 meV out of a total energy of 28 meV. This 
substantial anisotropy suggests that a single site term 
should be added to the Kitaev-Heisenberg Hamiltonian, 
probably resulting in a rather different phase diagram. 
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With all these caveats, it is still instructive to analyze 
where Na2ir03 is to be found in the parametric space of 
the nnKH model. We make the following assumptions: 
(i) that the atomic orbitals are fully localized and the 
appropriate basis is given by pure j e ff = 1/2 orbitals; (ii) 
that the only hoppings relevant for magnetic interactions 
are pd hoppings, so that the only oxygen assisted Ir-Ir 
hoppings are specific t 2g — t% g hoppings between unlike 
orbitals, as outlined in Refs. 12151 and the t 2g — e g hop- 
pings given in Ref. (4) and (iii) that the only processes 
contributing to magnetic interactions are those listed in 
Ref. H 

Indeed, the fact that the experimentally observed mag- 
netic order is zigzag suggests that either the Heisenberg 
terms are exceptionally long ranged (t he 3rd neighbor 
exchange is comparable to the 1st onef ^, or that the 
Kitaev term is strong and antiferromagneticP^l. The for- 
mer suggestion is seemingly in contradiction with the fact 
that the calculated 3rd neighbor hoppings are substan- 
tially smaller than the 1st neighbor ones. This makes 
it impossible to explain the large 3rd neighbor exchange 
integral in terms of superexchange. However, there is 
a possibility, suggested in Ref. that the Ir electrons 
are itinerant over individual hexagons, which makes mag- 
netic interactions naturally long ranged, and not directly 
related to the hopping integrals. 

The second suggestion, which is the one we will focus 
on in what follows, was proposed in Ref. HI namely that 
of an antiferromagnetic Kitaev term. If strong enough, 
this could explain the observed magnetic order. Below 
we consider the expressions presented in Ref. H] and sub- 
stitute the unknown variables with ab initio-dei'wed pa- 
rameters. 

Chaloupka et al™ discuss four relevant processes con- 
tributing to the exchange interactions in Na2lrC>3: (1) 
Direct hopping t\ a between nearest neighbor Ir t 2g or- 
bitals contributing with a term I± = (|iio-) 2 /U to the 
Heisenberg term, where U is the Coulomb repulsion be- 
tween t 2g electrons. 

(2) Interorbital nearest neighbor Ir t2 g -e g hopping via 
intermediate oxygens t\, with t% = tpdo-ipdw/A, where A 
is the charge-transfer energy (the difference between the 

A P f 

O p and Ir d levels) contributing with a term 1% = | -± ^ 
both to the Kitaev and Heisenberg terms, but with the 
opposite signs. Here U is the excitation energy associated 
with the t 2g -e g hopping i.e. it also includes crystal field 
splitting, U = U + lODq. Jh is the Hund's rule coupling 
between t 2g and e g electrons. 

(3) Oxygen-assisted hopping between two nearest 
neighbor Ir t 2g orbitals t\ o contributing with a term 

I 3 = | t jj 1 4& to the Kitaev term, where Jh is the Hund's 
rule coupling between t 2g electrons, and, we remind, 
*io=^ Ar /A. 

(4) Oxygen- 2p - Iridium-5d charge transfer contribut- 
ing with a term J 4 = ^ ^+^-3^ + 3( 2A+^-j ii ) + 
3(2A+u +2J ) — 3]' wnere U p and J p are, respectively, 



parameter 


value (eV) 


meaning 




0.03 


direct Ir-Ir hopping 


tio 


0.27 


O assisted Ir-Ir hopping 


t 


0.38 


Ir t 2g - e g hopping 


tdp 


0.57* 


Ir-0 hopping 


A p d 


2.4 


charge transfer energy 


Jh 


0.5 


Ir t2 g Hund's rule coupling 


Jh 


0.5 


Ir tig — e g Hund's rule coupling 



TABLE III: DFT calculated values of transfer integrals and 
charge transfer energy for Na2lrC>3 and estimates of Hund's 
rule coupling strength as described in the text. 

the Hubbard repulsion and the Hund's rule parameter 
for oxygen. This expression was derived by G. Khali- 
ullirP^and is worth some additional discussion. The first 
three terms correspond to processes where two holes of 
the same or of opposite spins meet at an oxygen atom. 
Neglecting J p , one gets simply — ^ 2 , which re- 
flects the fact that if the Ir atoms have opposite spins 
one can create an intermediate state with two holes on 
the same oxygen orbital, which lowers the total energy. 
The last term appears due to ring exchange, with an in- 
termediate state where two holes are located on different 
oxygens. This process is only allowed when the ground 
state is FM, and only if the ground state hole is in an a\ g 
or Jeff = 1/2 state, but not for pure t 2g orbitals. How- 
ever, contrary to a common misconception, J p is large, 
between 1.2 and 1.6 eV. We have estimated U p and J p , 
using the technique described in Ref. 1371 and obtained 
Up = 2.7 and J p = 1.6 eV, consistent with earlier DFT 
estimates^. For non-relativistic orbitals it is compara- 
tively straightforward to account for the Hund's rule cou- 
pling on O, but for relativistic orbitals it becomes more 
tedious. 

If we expand I4 in both U p and J p , then 1 4 w 

— ^o- U 2^2 P ■ This expression shows that the effect of U p 
and J p largely cancels and this term appears to be unim- 
portant (note though that if J p is entirely neglected, as 
in Ref. |H this proposition becomes more questionable). 

Summarizing the above terms into a single expression, 
Eq. ([!]) can be written as: 

H = (27 2 - h + 2/4) S?S] + (h -I 2 - h) Si ■ Sj. (17) 
2K J 

This model has a zigzag magnetic ground state^ if the 
Kitaev term is antiferromagnetic (AFM) and the Heisen- 
berg term is ferromagnetic (FM), with K > 0, J < and 
0.038 < \ J/K\ < 3.02. 

In Table |III| we provide the parameter values relevant 
for Na2lr03, as obtained from our DFT results. Note 
that the t\ parameter was assumed to be 2t\ o in Ref. El 
while in the calculations (DFT calculations are usually 
very reliable in this aspect) fi/iio 1S 1-4- However, using 
the ratio of 2 hardly changes any conclusions. 
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We present our results in Figures [9] and 10 In Fig. [9] 
we show the calculated values of K and J as a function 
of two variables: the x axis is the Hubbard U associ- 
ated with the upper Kramers doublet, and the y axis 
is the energy U, associated with exciting an individual 
electron from the upper ti g to an average e g state. The 
Hubbard U for 5d electrons is, generally speaking, 1.5 to 
2 eV. However, in this case it is additionally screened by 
the e g electrons, and also reduced by hybridization (cf. 
Na^CoOj^and Fe pnictides^) . Experimental estimates 
of the Hubbard U defined as the energy cost for excit- 
ing electrons across the insulating gap (which is th e def- 
inition relevant to superexchange) yield 0.3-0.5 e V n l 33 l 
Additionally, LDA+U calculations with U ~ 2 eV yield 
an excitation gap of the same order. We conclude that 
the realistic range of this parameter is 0.5-2 eV, with the 
smaller values more likely. 

For the second parameter, U, DFT calculations give ~ 
2.5 eV. This should be considered as a lower bound since 
DFT tends to slightly overestimate the orbital overlap 
and crystal fields, and misses the effects of the t 2g — e g 
Hubbard interaction. One can thus limit the physically 
admissible range in the region 2.5 eV < U < 3 eV. 

In Figure [10] we show the phase diagram in the space 
of the two parameters above. Several observations are in 
place: (1) While there is a zigzag phase in this diagram, 
it is very far removed from the range of the parameters 
that can be called physical, 0.5 eV < U < 2 eV, 2 eV< 
U < 3 eV (even though in the above estimate we have lib- 
erally stretched the admissible range in favor of a zigzag 
phase). In fact, the zigzag regime appears only when 
U < 0.6J7, i.e. when the Hubbard gap is larger than the 
e g — t2 g splitting, a rather unlikely proposition. (2) In 
the physical range of parameters, the ground state is ei- 
ther ferromagnetic or the spin liquid phase. It is rather 
curious that the very narrow slivers of the phase space in 
the J, K coordinated are transformed into a very large 
range in the U, U space. 

It is also worth mentioning that in order to explain 
the experimental data of Refs. l7l8l one needs not only the 
ground state to be zigzag, but also that K be several 
times larger than |J|; Chaloupka et aiPused K = 10.44 
and J = —4.01 eV. However, a closer look at the expres- 
sions in their work reveals that K + J = Ii — is/2, which 
does not depend on U and is always negative. Thus the 
two equalities above cannot be satisfied simultaneously 
for any choice of parameters, be they physical or not. 
Moreover, the values of J and K used in Ref. H] can only 
be obtained if U < 0.2 eV, which is clearly an impossible 
regime. 



B. Long-range exchange 

As mentioned above, an alternative interpretation of 
the experimental results, given in Refs. l7|8l is in terms of 
sizable 2nd and 3rd neighbor exchange constants, com- 
parable to the nearest neighbor exchange. In this picture 



the Kitaev term may or may not play a role, but this 
role is not decisive in establishing the observed magnetic 
order. Given that the calculated hopping amplitudes 
(Table I) are clearly dominated by the nearest neighbor 
terms, standard superexchange cannot explain such long 
range interactions. 

However, it is important to remember that in the op- 
posite, itinerant limit every electron is fully delocalized 
over a hexagon and, as such, is equally sensitive to the 
mean field magnetization pattern on the 1st, 2nd or 3rd 
nearest neighbors. As discussed in our earlier worlP, the 
zigzag order, as compared to the stripy one, results in a 
sizable pseudogap at the Fermi level even without a Hub- 
bard U. This creates an energy gain that cannot be cast 
in a form of nearest neighbor interaction, as it depends 
on the magnetization pattern over an entire hexagon. 

We are far from stating that the superexchange Hamil- 
tonian outlined in Ref. |4] is irrelevant and an itinerant 
description will give the final answer to all questions re- 
garding the magnetism in Na2lrC>3. However, relying 
solely on the localized picture and, correspondingly, on 
the nnKH model, is, apparently, inadequate. 



VII. CONCLUSIONS 

In summary, we have performed an extensive investiga- 
tion of the electronic properties of Na 2 Ir03 in the frame- 
work of non-relativistic and relativistic density functional 
theory calculations and derived by means of the Wan- 
nier function projector method, the corresponding mi- 
croscopic parameters. We resolved the following open 
questions: (1) By considering various idealized crystal 
structures for Na2lrC>3 we could disentangle the effect of 
each of the structural distortions present in this system 
and concluded that it is the joint effect of these distor- 
tions that constructively enhances the intrahexagon effec- 
tive hopping parameters and suppress the interhexagon 
ones favoring the formation of quasi-molecular orbitals. 
(2) We modelled the relativistic DFT results in terms of 
a tight-binding model including the spin-orbit coupling 
term and analyzed the electronic properties of Na2lr03 
in terms of two complementary descriptions, the (itiner- 
ant) quasi-molecular basis and the (localized) relativistic 
j e s basis. We observed that the behavior of Na2ir03 
lies in between the fully itinerant and the fully localized 
description and that a quasi-molecular orbital descrip- 
tion keeps its character even at large values of the spin- 
orbit coupling strength. (3) We showed that XAS and 
RIXS observations can be well understood within an itin- 
erant description of Na2lr03. (4) Finally, we provided ab 
initio-derived estimates for the parameters appearing in 
the Kitaev and Heisenberg terms in Na2lrC>3 and found 
that the recently proposed nnKH model (see Section VI I , 
even though it is a very interesting model per se, is un- 
fortunately not realistic for Na2lr03. 
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(a) K(meV) (b) J(meV) 




U (eV) U (eV) 

FIG. 9: (Color online) Variation of (a) Kitaev parameter K and (b) Heisenberg exchange coupling J with onsite Coulomb 
coupling strength U and Ir t2 g -e g excitation energy U. Positive values refer to antiferromagnetic, negative to ferromagnetic 
values of K and J. The other parameters entering the K and J are given in Table [TTT| 



FIG. 10: (Color online) Phase diagram of the Kitaev- 
Heisenberg model for Na2lr03 with parameters determined 
following Ref. [4] The calculated exchange integrals are func- 
tions of the Mott-Hubbard gap U and the cubic crystal field 
splitting U. The contours mark isolines of the ratio K/J. 
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VIII. APPENDIX. IDEALIZED Na 2 Ir0 3 
CRYSTAL STRUCTURES AS USED IN WIEN2k 



A. Experimental structure from Ref. 8 



Space group C2/m (No. 12) 
a = 5.4269 A, b = 6.4104 A, c = 9.3949 A, 7 = 124.12° 
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0.0 


0.0 


0.1670 


01 


0.4590 


0.2110 


0.1780 


02 


0.0070 


0.7960 


0.0 



B. Idealized structure 5*1 



Space group C2/m (No. 12) 
a = 5.4501 A, b = 6.4411 A, c = 9.4399 A, 7 = 125.26° 



Atom x y z 

Nal O0 0.0 05~ 

Na2 0.5 0.5 0.0 

Na3 0.5 0.5 0.3333 

Ir 0.0 0.0 0.1667 

01 0.4646 0.2097 0.1785 

02 0.0000 0.7903 0.0 



Atom x y z 
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C. Idealized structure S2 



Space group C2/m (No. 12) 
a = 5.4501 A, b = 6.4190 A, c = 9.4399 A, 7 = 124.47° 



Atom x y z 

~NaI (10 O0 0lT~ 

Na2 0.5 0.5 0.0 

Na3 0.5 0.5 0.3333 

Ir 0.0 0.0 0.1667 

01 0.4606 0.1909 0.1667 

02 0.0000 0.8091 0.0 



Note that due to the necessity of using a monoclinic 
angle 7 in WIEN2k, the Ir honeycomb layers in the 
Na2lr03 unit cells presented above are parallel to the 
ac plane. Accordingly, within this convention the vector 
of the Bloch factors in Eqs. Q and ^ is given by 



D. Idealized structure S3 



T: 



M=1...6 



(k) 



(l,e- ifc * s ,e ifc * 5 ,e i(fc '' 



-k x a) ik z c - 



Space group C2/m (No. 12) 
a = 5.0658 A, b = 5.9869 A, c = 8.7743 A, 7 = 125.26° 
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1 D. J. Kin 


0.0000 
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0.0 

i, T. Morita, II. 



), 

where a and c are the lengths of the two primitive lattice 
vectors lying in the ac plane. Here, one explicitly ac- 
counts for the choice of WIEN2k of the actual positions 
of the two Ir atoms in the primitive unit cell, which are, 
e. g., (-0.167,0,0.167) and (-0.833,0,0.833) in the experi- 
mental Na2lr03 structure. 



Takagi, T. Arima, Science 323, 1329 (2009). 

J. Chaloupka, G. Jackeli, and G. Khaliullin, Phys. Rev. 

Lett. 105, 027204 (2010). 

J. Reuther, R. Thomale, and S. Trebst, Phys. Rev. B 84, 
100406(R) (2011). 

J. Chaloupka, G. Jackeli, and G. Khaliullin, 



arXiv:1209.5100 (unpublished). 



I. I. Mazin, H. O. Jeschke, K. Foyevtsova, R. Valentf, and 
D. I. Khomskii, Phys. Rev. Lett. 109, 197201 (2012). 

6 The spin-orbit splitting is given by the difference of the 
eigenvalues A/2 an — A in the j c g space. 

7 Y. Singh, S. Manni, J. Reuther, T. Berlijn, R. Thomale, 
W. Ku, S. Trebst, P. Gegenwart, Phys. Rev. Lett. 108, 
127203 (2012). 

8 S. K. Choi, R. Coldea, A. N. Kolmogorov, T. Lancaster, 
I. I. Mazin, S. J. Blundell, P. G. Radaelli, Y. Singh, P. 
Gegenwart, K. R. Choi, S.-W. Cheong, P. J. Baker, C. 
Stock, and J. Taylor, Phys. Rev. Lett. 108, 127204 (2012). 

9 From these considerations, one can think of other materials 
with basically the same structure, in which 1/3 of the in- 
plane foreign atoms M in Na(M!/3lr 2 /3)02 can be any 
metal M. Such a program was recently realized by the 
group of R. J. Cava (privat e communication). 

10 Y. Singh and P. Gegenwart, Phys. Rev. B 82, 064412 
(2010). 

11 R. Comin, G. Levy, B. Ludbrook, Z.-H. Zhu, C. N. Veen- 
stra, J. A. Rosen, Y. Singh, P. Gegenwart, D. Strieker, J. 
N. Hancock, D. van der Marel, I. S. Elfimov, A. Damascelli, 
Phys. Rev. Lett. 109, 266406 (2012). 

12 P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvas- 
nicka, and J. Luitz 2001 WIEN2k, An Augmented 
PlaneWave+LocalOrbitals Program for Calculating Crystal 
Properties (Karlheinz Schwarz, Techn. Universitat Wien, 



Austria) . 

13 J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett. 
77 3865 (1996). 

14 S. Cottenier, Density Functional Theory and the fam- 
ily of (L)APW-methods: a step-by-step introduction 
(Instituut voor Kern- en Stralingsfysica, K.U.Leuven, 
Belgium), 2002, ISBN 90-807215-1-4 (to be found at 
http : //www . wien2k. at /reg_user /textbooks). 

15 M. Aichhorn, L. Pourovskii, V. Vildosola, M. Ferrero, O. 
Parcollet, T. Miyake, A. Georges, S. Biermann, Phys. Rev. 
B, 80, 085101 (2009). 

16 J. Ferber, K. Foyevtsova, H. O. Jeschke, R. Valentf, 
|arXiv:1209.4466| (unpublished) . 

17 The description of the orbital character in the present work 
is done according to the local cubic coordinate system if 
not explicitly specified otherwise. 

18 K. Ko epernik and H. Eschri g, Phys. Rev. B 59. 1743 
(1999) ; Jhttp: / /www.FPLO.de| 

19 In Fig. [6] we visualize real- valued combinations of degener- 
ate E 2u or Ei g QMOs. 

20 We compared our prepared idealized structures to those 
provided by Radu Coldea (private communication) and the 
agreement between both sets of structures is very good. 

21 C. Martins, M. Aichhorn, L. Vaugier, S. Biermann, Phys. 
Rev. Lett. 107, 266404 (2011). 

22 A. Shitade, H. Katsura, J. Kunes, X.-L. Qi, S.-C. Zhang, 
and N. Nagaosa, Phys. Rev. Lett. 102, 256403 (2009). 

23 Note that these states are eigenstates of Eq. (jfij) after diago- 
nalizing the matrix Eq. (JsJ) and performing an appropriate 
rotation. 

24 Orthorhombic distortion lifts the degeneracy of the E\ g 
and E2u states combining them into real functions. 

25 In the ideal case, one would have, e. g., B\ u f= 
(-0.408, -0.408, -0.408, 0.408, 0.408, 0.408, 0, 0, 0, 0, 0, 0). 



17 



D. Haskel, G. Fabbris, M. Zhernenkov, P. P. Kong, C. Q. 
Jin, G. Cao, and M. van Veenendaal, Phys. Rev. Lett. 109, 
027204 (2012). 

L. C. Chapon and S. W. Lovesey, J. Phys.: Condens. Mat- 
ter 23, 252201 (2011). 

S. W. Lovesey and A. N. Dobrynin, J. Phys.: Condens. 
Matter 24, 328201 (2012). 

S. W. Lovesey, D. D. Khalyavin, P. Manuel, L. C. Chapon, 
G. Cao and T. F. Qi, J. Phys.: Condens. Matter 24, 496003 
(2012). 

G. van der Laan and B. T. Thole, Phys. Rev. Lett. 60, 
1977 (1988). 

B. T. Thole, and G. van der Laan, Phys. Rev. B 38 3158 
(1988). 

J. P. Clancy, N. Chen, C. Y. Kim, W. F. Chen, K. W. 
Plumb, B. C. Jeon, T. W. Noh, and Y.-J. Kim, Phys. Rev. 



B 86, 195131 (2012). 

H. Gretarsson, J. P. Clancy, X. Liu, J. P. Hill, E. Bozin, 
Y. Singh, S. Manni, P. Gegenwart, J. Kim, A. H. Said, 
D. Casa, T. Gog, M. H. Upton, H.-S. Kim, J. Yu, V. M. 
Katukuri, L. Hozoi, J. van den Brink, and Y.-J. Kim, Phys. 
Rev. Lett. 110, 076402 (2013). 

Y. Yu, L. Liang, Q. Niu, S. Qin, 87, 041107(R) (2013). 

G. Khaliullin, private communication. 

I. I. Mazin and D. J. Singh, Phys. Rev. B 56, 2556 (1997). 
A. G. Petukhov, I. I. Mazin, L. Chioncel and A. I. Lichten- 
stein, Phys. Rev. B 67, 153106 (2003). 

H. Ishida, M. D. Johannes, and A. Liebsch, Phys. Rev. 
Lett. 94, 196401 (2005). 

See for instance, M. Aichhorn, L. Pourovskii, and A. 
Georges, Phys. Rev. B 84, 054529 (2011). 



